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We report on v e and v e appearance in and beams using the full MINOS data sample. The 
comparison of these v £ and i/ e appearance data at a 735 km baseline with On measurements by 
reactor experiments probes 8, the 623 octant degeneracy, and the mass hierarchy. This analysis is 
the first use of this technique and includes the first accelerator long-baseline search for P M — > P e . Our 
data disfavor 31% (5%) of the three-parameter space defined by S, the octant of the 823, and the 
mass hierarchy at the 68% (90%) C.L. We measure a value of 2sin 2 (2(9i3)sin 2 (#23) that is consistent 
with reactor experiments. 
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PACS numbers: 14.60.Pq, 14.60.Lm, 29.27.-a 

The neutrino oscillation phenomenon is successfully 
modeled by a theory of massive neutrino eigenstates that 
are different from the neutrino flavor eigenstates. These 
sets of eigenstates are related by the PMNS matrix [T] 
which is commonly parameterized by three angles, 9ij, 
and a CP-violating phase, S. 

The values of 812 and 823 have been measured [2H!] 
with indications that 6*23 is not maximal 5 7. The final 
angle, 813, is now known to have a nonzero value from 
measurements by reactor experiments [HrllOj. the mea- 
surement by the T2K accelerator experiment, and 
from earlier MINOS results [121113]. 

Despite these accomplishments, the value of S is still 
unknown, as is the ordering of the neutrino masses, which 
is referred to as the neutrino mass hierarchy. Much of 
the attention in the neutrino community is now focused 
on resolving these unknowns. The mass hierarchy is not 
only a fundamental property of neutrinos but also has a 
direct impact on the ability of ncutrinolcss double beta 
decay searches to state definitively whether the neutrino 
is its own antiparticle [14] . Reactor experiments make 
a pure measurement of #13, whereas the — > v e and 
— > v e appearance probabilities measured by acceler- 
ator experiments such as MINOS depend on the value 
of S and sin 2 (6>23). In addition, the long-baseline of MI- 
NOS means that interactions between neutrinos and the 
matter of the Earth make the appearance probabilities 
dependent on the neutrino mass hierarchy jTS], [TB] . 

We report the result from the search for v e (p e ) appear- 
ance in a (Dp) beam using the full MINOS data sample. 
This result uses an exposure of 10.6 x 10 20 protons-on- 
target taken with a v beam and an exposure of 3.3 x 10 20 
protons-on-target taken with a v beam. The neutrino 
sample is 30% larger than the sample used for the pre- 
vious MINOS results on this topic [T3]. This analysis 
represents the first long-baseline search for i> M — > v e ap- 
pearance and places new constraints on #13 and on a 
combination of d, 823, and the neutrino mass hierarchy. 

In the MINOS experiment [T7], neutrino oscillation is 
studied with the NuMI beamline [18 by measuring neu- 
trino interactions in two detectors. The Near Detector 
(ND), which has a fiducial mass of 29 tons, is at a dis- 
tance of 1.04 km from the production target and is used 
to determine the composition of the beam before the neu- 
trinos have oscillated. The Far Detector (FD) , which has 
a fiducial mass of 3.8 kilotons, is at a distance of 735 km 
from the production target and is used to measure the 
change in the neutrino flavor composition of the beam. 
In both the v and v beam modes, the NuMI beam has an 
energy spectrum that is peaked at 3 GeV. At the ND, the 
neutrino flavor composition of the neutrino interactions, 
as determined by a combination of simulation and mea- 
surement, is found to be 91.7% 7.0% and 1.3% 



v e and v e for the v beam mode and 58.1% v^, 39.9% v^, 
and 2.0% v e and i/ e for the v beam mode. 

Both detectors are magnetized tracking calorimeters 
consisting of alternating planes of 2.54 cm thick steel 
and 1 cm thick scintillating plastic [T7]- The scintilla- 
tor planes are segmented into 4.1 cm wide strips with 
wavelength-shifting fibers embedded in the strips to col- 
lect light for readout by multi-anode photomultiplier 
tubes. 

In the MINOS data sample, the flavor of a neutrino is 
determined only for charged-current (CC) interactions. 
Vp-QG and D^-GG interactions are identified by the pres- 
ence of a long muon track that extends beyond a cluster 
of energy depositions that are consistent with hadronic 
activity at the interaction vertex. Neutral-current (NC) 
interactions are identified by the energy depositions asso- 
ciated with hadronic activity. v e -CC and P e -CC interac- 
tions produce an electromagnetic shower that typically 
leaves a compact cluster within 6 to 12 planes. This 
analysis does not distinguish between v e -CC and v e -CC 
interactions. 

The sample of events classified as ^ e -CC and v e -CC 
interactions contains a background of interactions with 
similar topology as required for v e -GG and v e -CC classi- 
fication. NC interactions with a significant electromag- 
netic component and v^-GC or v^-CG interactions in 
which the muon track is not easily identified make up 
the majority of the background. Smaller contributions 
to the background arise from v T -CQ, and v r -GG interac- 
tions. In addition to backgrounds that mimic v e -CC and 
P e -CC event topologies, intrinsic v e and D e components 
of the NuMI beam must be taken into account. 

Candidate ^ e -CC and v e -CC events are required to fall 
within a fiducial volume and to be coincident in time and 
direction with the NuMI beam. We require the events 
to have shower-like topologies by rejecting events with 
tracks that are longer than 25 planes or extend more 
than 15 planes from a shower edge. In addition, recon- 
structed events must have at least five consecutive planes 
with deposited energy above a threshold; this threshold 
is defined as half of the energy deposited by a minimum 
ionizing particle. We require the events to have a recon- 
structed energy between 1 and 8 GeV where most of the 
v e and T> e appearance is expected. 

We further classify the events in this pre-selected 
sample of shower-like events by using a library- event- 
matching (LEM) algorithm PH [20]. Within the LEM 
algorithm, the topology of energy depositions that char- 
acterize the event is compared to a library of simulated 
signal and background events. Separate libraries are used 
for the events in the v beam mode and T> beam mode. The 
50 best-matching events in the library are collected and 
used to produce three variables. These variables are the 
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FIG. 1: Distributions of ollem- The plots in the left column 
correspond to the v beam mode. The plots in the right col- 
umn correspond to the v beam mode. The top row shows the 
distributions for ND selected events with a band about the 
simulation representing the systematic uncertainty. The bot- 
tom row shows the distributions for the predicted FD back- 
ground and signal multiplied by 10 with 2sm 2 (2#i3)sin 2 (6 l 23) 
= 0.1, 5 = 0, and a normal mass hierarchy. 



fraction of best-matching library events that are u e -CC 
or £/ e -CC, the average inelasticity of the best-matching 
v e -CC or D e -CC library events, and the average frac- 
tion of the energy depositions that overlap between the 
test event and the best-matching v e -CC or D e -CC library 
events. These three variables and the reconstructed neu- 
trino energy of the test event are then used as an input 
into an artificial neural network. The output value from 
the neural network is used to discriminate between signal 
and background events. This discriminant variable is re- 
ferred to as a lem an d is shown in Fig. [T] Signal events 
have a value near one, while background events cluster 
near zero. The maximum sensitivity to v e and v e appear- 
ance is obtained by analyzing events with ollem > 0.6. 

Following the selection of v e -CC and D e -CC candidate 
events, the ND data are used to study the rate of back- 
ground from NC, z/^-CC and P^-CC and intrinsic beam 
v e -CC and v e -GG interactions. The NuMI beam can be 
tuned to produce different energy spectra. Among these 
different beam configurations, the relative contributions 
of the various backgrounds change in a well-understood 
way. By measuring the total of the three backgrounds in 
three different beam configurations, the relative amounts 
of the individual backgrounds can be deduced [21] . 

We use the measurement of the ND backgrounds to de- 
rive the FD background predictions for the data samples 
in the v beam mode and in the v beam mode. For each 
sample, we divide simulated FD events into bins of energy 



Systematic Effect Uncertainty Uncertainty 
v mode V mode 


Energy Scale 


2.7% 


3.0% 


Normalization 


1.9% 


1.9% 


v r cross-section 


1.7% 


2.0% 


All Others 


0.8% 


2.5% 


Total Systematic 


3.8% 


4.8% 


Total Statistical 


8.8% 


23.9% 



TABLE I: Systematic uncertainty on the FD background pre- 
diction for events with a value of ollem > 0.6. Effects listed 
under "All Others" include the neutrino flux, cross sections, 
detector modeling, and background decomposition. 



and a lem and correct each FD background component, 
bin-by-bin, by multiplying it by the measured ND ratio 
of data to simulated events for that background. Since 
the ND data sample does not contain v T -CC and D T -CC 
events from oscillation, we estimate the FD contribution 
from this small background component through simula- 
tion and a correction based on the observed ND ^ M -CC 
and v^-CC spectra. 

The sources of systematic uncertainty that affect the 
background prediction are given in Table [TJ The effect 
of each source of uncertainty is evaluated by producing 
simulated ND and FD event samples that are modified 
according to the estimated size of each systematic effect. 
These modified samples are used to produce an altered 
FD background prediction for the systematic effect in 
question. We take the resulting difference between the 
nominal and modified prediction as the systematic un- 
certainty on the background prediction. The systematic 
effect that results in the largest reduction in sensitivity is 
a 2.0% uncertainty on the relative energy scale between 
the ND and FD. 

With the absence of a i/ e -CC and D e -CC signal in 
the ND, the signal selection efficiency cannot be ex- 
trapolated from the ND events in the same way as the 
background estimate. Therefore, to evaluate the signal 
efficiency, we select a sample of well identified v^-CC 
events [221 US] > remove the energy depositions that are 
associated with the muon track |24| , and insert the simu- 
lated energy depositions of an electron with an identical 
three-momentum |25j . This method effectively turns a 
well identified sample of z^-CC and D^-CC data events 
into a sample of v e -CC and P e -CC data events. For the 
v beam mode (y beam mode) data sample, we find the 
expected number of FD signal events with ollem > 0.6 
and the associated systematic uncertainty to be 33.7±1.9 
(3.9 ± 0.2) assuming sin 2 (26»i 3 ) = 0.1, (5 = 0, 6» 23 = tt/4, 
and a normal mass hierarchy. This corresponds to an 
identification efficiency of (57.4 ± 2.8)% for the v beam 
mode and of (63.3 ± 3.1)% for the v beam mode. The 
systematic uncertainties are evaluated in a way that is 
similar to the evaluation of the background systematics 
by using simulated samples that have been altered by a 
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Event Type 


^ beam 


v beam 




mode 


mode 


NC 


89.4 


13.9 


f M -CC and P^-CC 


21.6 


1.0 


Intrinsic i/ e -CC and P e -CC 


11.9 


1.8 


v T -\^K_j ana y T -UL 


4.0 


U.6 


I/jj — S> fe-CC 


33.0 


0.7 


— > i/ e -CC 


0.7 


3.2 


Total 


161.4 


21.4 


Data 


152 


20 



MINOS Far Detector Data 



TABLE II: Expected FD event yields for events with a value 
of aLEM > 0.6, assuming sin 2 (20i3) = 0.1, 5 = 0, 023 = 7i"/4, 
and a normal mass hierarchy. 



systematic effect. 

Events with aLEM < 0.5 are insensitive to v e and v e 
appearance. These events are therefore used in a separate 
study to validate the analysis procedure. ND events with 
&lem < 0.5 are used to predict FD event yields, which 
are found to agree with the FD data to within 0.3 (0.6) 
standard deviations of the statistical uncertainty for the 
data sample in the v (y) beam mode. 

Events with ulem > 0.6 are selected for further anal- 
ysis in the v beam mode and in the v beam mode. The 
expected and observed event counts in these samples are 
shown in Table |Tij The observed FD reconstructed energy 
spectra, in bins of olem, are shown for the candidate 
events in Fig. [2j Assuming a three-flavor neutrino oscil- 
lation probability that includes matter effects [IB], we 
simultaneously fit the data from the v beam mode and v 
beam mode samples for the value of 2sin 2 (20i 3 )sin 2 (0 23 ) 
while the value of the mass hierarchy and 5 are held 
fixed. The fit is performed using the 15 bins formed 
by three bins of chlem and five bins of energy. This 
procedure is performed for all values of 5 and both 
mass hierarchies, and the resulting confidence intervals, 
calculated using the Feldman-Cousins technique [2"B] . 
are shown in Fig. [3j The values of the oscillation 
parameters used in the fit are taken from previous mea- 
surements dS] and are set to sin 2 (20 23 ) = 0.957 tomh 
|Am§ 2 | = (2.39 to°l) x 10~ 3 eV 2 , 9 12 = 0.60 ± 0.02, 
and Amfi = (7.59 x 10~ 5 eV 2 . The full set of 

statistical and systematic uncertainties on the prediction 
are taken into account when constructing the contours. 

Assuming a normal mass hierarchy, S = 0, and 
#23 < tt/4, we find that the data allow for values 
of 0.01 < 2sin 2 (20i 3 )sin 2 (0 23 ) < 0.12 at 90% C.L. with 
the best-fit value of 2sin 2 (20 13 )sm 2 (0 23 ) = 0.051+^. 
Assuming an inverted mass hierarchy, 6 = 0, and 
#23 < 7r/4, we find that the data allow for values of 



0.03 < 2sin 2 (26>i 3 )sin 2 (0 23 ) < 0.18 at 90% C.L. with the 
best-fit value of 2sin 2 (26>i 3 )sin 2 (0 23 ) = 0.093^;^. The 
best-fit values show very weak dependence on the choice 
of octant for 6* 23 . 

We are further able to place constraints on the 
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Reconstructed Energy (GeV) 

FIG. 2: The reconstructed energy distributions for three 
OLhEM ranges. The events with energy greater than 5 GeV 
are combined into a single bin for the fits. The vertical bars 
through the data points denote statistical uncertainties. The 
signal predictions assume sin (20i 3 ) = 0.051, Am| 2 > 0, 8 = 
0, and 023 = 7r/4. The plots in the left column correspond to 
data collected in the v beam mode. The plots in the right 
column correspond to data collected in the v beam mode. 



value of <5, the octant of 23 , and the neutrino mass 
hierarchy by incorporating the current knowledge of 
sin 2 (2#i 3 ) = 0.098 ± 0.013 that we calculate from recent 
reactor data [5HTU]. Fig. [4] shows the likelihood for our 
data as a function of 6 for the four possible combinations 
of mass hierarchy and the octant of # 23 . The full set 
of statistical and systematic uncertainties on the predic- 
tion are taken into account when calculating the likeli- 
hood, as are the uncertainties on the oscillation param- 
eters. This analysis represents the first result by a long- 
baseline experiment to use a combination of — > v e 
and — > v e appearance data, with external reactor con- 
straints on 0i 3 , to probe (5, the #23 octant degeneracy, and 
the mass hierarchy. Assuming # 23 > 7r/4 (# 23 < 7r/4), the 
data prefer an inverted hierarchy at 0.63 (0.04) units of 
— 2AlnL. Furthermore, as is indicated by the horizontal 
lines in Fig. [4j our data disfavor 31% (5%) of the three- 
parameter space defined by 5, the octant of the 23 , and 
the mass hierarchy at the 68% (90%) C.L. This analysis 
demonstrates the potential discriminating power that can 
be achieved with the combination of reactor and — > v e 
and — > v e appearance data. 

In conclusion, we have presented the results of v & and 
v e appearance in and beams from the full MINOS 
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FIG. 3: The 68% and 90% confidence intervals of allowed 
values for 2sin 2 (2#i3)sin 2 (#23) as a function of 5 for the two 
mass hierarchies. 



data sample. We have used these data to place new con- 
straints on the mixing angle #13 and have demonstrated 
how such data will be used in the future to break the 
degeneracy in the appearance probability created by the 
ambiguity in the octant of 023, the neutrino mass hierar- 
chy, and the value of the CP-violating phase S. 
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FIG. 4: The resulting values of the likelihood L, shown here 
as — 2AlnI/, from a fit of 8 to our data using constraints from 
reactor experiments [8til0|. assuming various values of the 
mass hierarchy and the sign of #23 — 7r/4. The difference is 
taken with respect to the best-fitting solution. Values above 
the horizontal dashed lines are disfavored at either 68% or 
90% C.L. 
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